Lissajous curves and semiclassical theory: The two-dimensional harmonic oscillator 
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The semiclassical treatment of the two-dimensional harmonic oscillator provides an instructive 
example of the relation between classical motion and the quantum mechanical energy spectrum. 
We extend previous work on the anisotropic oscillator with incommensurate frequencies and the 
isotropic oscillator to the case with commensurate frequencies for which the Lissajous curves appear 
as classical periodic orbits. Because of the three different scenarios depending on the ratio of its 
frequencies, the two-dimensional harmonic oscillator offers a unique way to explicitly analyze the 
role of symmetries in classical and quantum mechanics. 
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I. INTRODUCTION 

In 1940 J. M. Jauch and E. L. Hill discussed the re- 
lation between symmetries and degeneracies in the spec- 
tra of quantum systems. As examples they considered 
the hydrogen problem, the two-dimensional Kepler prob- 
lem, and the isotropic harmonic oscillator, but the two- 
dimensional anisotropic harmonic oscillator resisted their 
analysis. The authors concluded "We feel, however, that 
the difficulties encountered in this relatively simple prob- 
lem throw question on the true interpretation of classical 
multiply-periodic motions in quantum mechanics."— 

Although the problem is now understood on a group 
theoretical level j2i2i4 the remark by Jauch and Hill mo- 
tivates us to study the relation between the motion of 
a classical anisotropic oscillator with commensurate fre- 
quencies and its quantum mechanical energy spectrum. 
This link is provided by semiclassical theory^ which re- 
lates the density of states to properties of the classical 
periodic orbits. We will restrict the discussion to the 
two-dimensional case where harmonic potentials already 
provide three different scenarios for periodic orbits. 

For the isotropic harmonic oscillator all orbits are peri- 
odic and constitute a continuous family of which two rep- 
resentatives are depicted in Fig. QJa). Any two distinct 
periodic orbits can be deformed continuously into each 
other via other periodic orbits. The situation is different 
for an anisotropic harmonic oscillator with incommensu- 
rate frequencies where only the two normal modes shown 
in Fig. QJb) correspond to periodic orbits. The most in- 
teresting classical dynamics is exhibited by the harmonic 
oscillator with commensurate frequencies of which an ex- 
ample is shown in Fig. [ljc). Here, Lissajous curves ap- 
pear as additional periodic orbits that form a continuous 
family. These trajectories, which will play an important 
role in the following, take their name after Jules An- 
toine Lissajous who gave an extensive discussion of these 
curves in a study aimed at a comparison of the relative 
frequencies of two tuning forks.— Previous work address- 
ing these curves should not be forgotten.— 

It has been proven that semiclassical theory gives the 
exact expression for the density of states of an isotropic 8 
as well as for the incommensurate harmonic oscillator. 9 
An interesting property of the harmonic oscillator is that 
the appropriate limit from incommensurate to commen- 
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FIG. 1: Examples of the periodic orbits of a two-dimensional 
harmonic oscillator. The dashed line corresponds to an 
equipotential line, (a) Isotropic harmonic oscillator, (b) 
Anisotropic harmonic oscillator with incommensurate fre- 
quencies. In this case only the two orbits shown are peri- 
odic, (c) Anisotropic harmonic oscillator with commensurate 
frequency ratio 1:2. In contrast to the case with incommen- 
surate frequencies, Lissajous curves become possible. 



surate frequencies yields the correct density of states even 
for the commensurate harmonic oscillator^ although 
only two periodic orbits are explicitly taken into account. 
Although taking such a limit considerably simplifies the 
calculation, it is unsatisfying that a full analysis of the 
density of states based on all the classical periodic or- 
bits is lacking. Furthermore, a complete account for the 
symmetries of the commensurate anisotropic harmonic 
oscillator is not possible without the Lissajous curves. 
As we will show, the proper treatment of all the peri- 
odic orbits of the commensurate anisotropic harmonic 
oscillator yields the exact density of states and the exact 
symmetry reduced densities of states. Our analysis thus 
closes the last lacuna left in the semiclassical treatment 
of the two-dimensional harmonic oscillator. This system 
now lends itself well to illustrating the relation between 
classical dynamics and the symmetry properties of the 
quantum mechanical energy spectrum. 



II. THE ENERGY SPECTRUM IN TERMS OF 
CLASSICAL MECHANICS: THE 
INCOMMENSURATE HARMONIC OSCILLATOR 



Among the three types of two-dimensional harmonic 
oscillators depicted in Fig. [IJ the relation between the 
quantum mechanical density of states and classical prop- 
erties is most readily established for the incommensurate 
harmonic oscillator. The Hamiltonian of a particle of 
mass m moving in a two-dimensional harmonic potential 
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with eigenfrequencies u\ and L02 is given by 

Pi m ? ? V\ m 1 2 /-.x 

We first use the fact that the total Hamiltonian H 
can be decomposed into two contributions Hi(qi,pi) 
and -^2(92,^2) describing two uncoupled one-dimensional 
harmonic oscillators. From the classical density of states 
of the ith normal mode, 



with 



Pi(E) 



2nn \ 2m 2 ) nuji 



(2) 



the classical density of states of the two-dimensional os- 
cillator is obtained by the convolution 



p(E)=[ dE'p l {E-E')p 2 {E') = 
Jo 



E 



(3) 



H 2 U)iUJ2 

The quantum mechanical energy spectrum is given by 
E V1<V2 =tiui(vi + ^) +^2(^2 + ^), (4) 

where ^1,^2 = 0,1,2,... The corresponding quantum 
mechanical density of states is 



V\ ,J/S=0 



(5) 



Equation ([5]) can be interpreted in terms of the classical 
properties of the incommensurate harmonic oscillator if 
it is appropriately re-expressed. We perform a Laplace 
transform of Eq. ([5]) and find the partition function of 
the incommensurate harmonic oscillator 
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Z p = dE p(E)exp(-fiE) 
Jo 



(6) 



4sinh y — - — J sinh ^ 



The inversion of the Laplace transform by means of a con- 
tour integration yields three kinds of contributions from 
the various poles on the imaginary axis. The double pole 
at (3 = gives rise to the classical density of states ([3]), 
and the poles at (3 = ±i2Trn/huji and f3 = ±i2Tm/huj2 
with n — 1,2,... lead to contributions responsible for 
the difference between the classical and quantum density 
of states. The latter can thus be expressed as 



p(E)=p(E)+S Pl (E)+5p2(E). 
After evaluating the residues, we obtain 

/ 2ttE tt 



(7) 



Spi(E) 



cos I 



2| sin(7rrw2/^i)| 



(8) 



oy 1 = 2r + 1 



W2 



(9) 



The quantities Sp2(E) and a r ^2 are expressed in the same 
way with the indices 1 and 2 interchanged. We have 
written the expressions in a form that anticipates the 
following discussion. In particular, we have introduced 
the periods Ti — 2n/uji of the two normal modes and 
phases a r .i- In Eq. (J9j) [x] denotes the integer part of x. 

We expect that it should be possible to relate the terms 
in 8 pi to the properties of the two periodic orbits shown 
in Fig. Hlb). As a motivation we remark that within the 
Feynman path integral approach 10 all paths joining two 
points in a given time contribute to the corresponding 
quantum propagator a phase factor exp(iS/h), where S 
is Hamilton's principal function of the path. For S 3> h 1 
we can apply the semiclassical approximation^^ and find 
that the path integral is dominated by the stationary 
points of S, that is, the classical paths. Within this 
approximation even the contribution of quantum fluc- 
tuations can be expressed in terms of properties of the 
classical motioni 12 ! 13 

The propagator can be related to the density of states 
by means of the Green function G(E) = (E — the 
Fourier transform of the propagator. From the relation 



p{E) 



- lim Im(Tr[G(£; + ie)l) 

7T e^0+ 
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7T e^0+ 



E 

Vl,V2 



(10) 
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we see that p{E) involves the trace over the Green func- 
tion, so that the initial and final points of the con- 
tributing classical paths coincide. Equation (flU|) con- 
tains contributions from very short orbits which yield 
p(E) and periodic orbits of finite length contribute to 
8p(E) = 8pi(E) + Sp2(E). The latter relation is given by 
the Gutzwiller trace formul a 12 ! 14 of which Eq. is an 
example. 

These considerations indicate that within the semiclas- 
sical approximation the density of states can be expressed 
in terms of classical properties, but such an expression 
cannot be expected in general. However, it was shown 
by Brack and Jain 9 that the incommensurate harmonic 
oscillator is one of the rare cases where the Gutzwiller 
formula yields the exact density of states. As a con- 
sequence, it should be possible to interpret the various 
terms in Eq. ([5]) by means of a classical periodic orbit. 
We will check now that this is the case. 

We have already mentioned that the prefactor T\ in 
Eq. ([8]) represents the period of the first normal mode. 
The argument of the cosine function in the numerator 
contains the classical action 

W 1 (E) = ldq lVl = 7 ^- (11) 
J Ux 

of the fundamental periodic orbit along the q\ axis for a 
given energy Hi = E; we sum over all repetitions counted 
by r. 



3 



The denominator is determined by the lateral stabil- 
ity of the periodic orbit as the following argument shows. 
The equations of motion of a harmonic oscillator can eas- 
ily be solved, and we find 



Pi(t) 



X<(t) 



with the matrizant 



Xi(i) 



«i(0) 



sin(wit) ' 
mtUi 

-vnwi sin(wjt) cos(w;i) 



(12) 



(13) 



We again consider the periodic orbit along the qi axis 
which is specified by the initial conditions 92(0) = and 
P2(0) = 0. We are interested in the evolution of a small 
lateral perturbation during r oscillations. After the time 
rXi, we find a deviation 



J*] = [x 2 i/-r 1 )-Ji] 



92(0) 
MO) 



(14) 



in the direction orthogonal to the periodic orbit with 
respect to the perturbed initial condition. The matrix 
X2(rTi) thus measures the stability during r repetitions 
of the first normal mode against a small initial perturba- 
tion. It enters the denominator of Eq. ([5]) by means of 
its property 



x/|det[X 2 (r-Ti) - i]j = 2| sin(7rrwaM)|. 



(15) 



The last quantity requiring interpretation is the 
Maslov index a rt i defined by Eq. ([9]). It collects all phases 
appearing in the evaluation of the stationary phase inte- 
gral of the first periodic orbit. Because the derivation of 
Eq. {9]) within a semiclassical approach is more demand- 
ing, we refer the reader interested in the details to the 
literature £ 

We started from an exact expression for the density 
of states of the incommensurate harmonic oscillator and 
arrived at a semiclassical interpretation. In contrast, in 
the following two sections, we will go the opposite way 
and first derive a semiclassical expression for the den- 
sity of states of the isotropic harmonic oscillator and the 
commensurate anisotropic harmonic oscillator by exploit- 
ing the system's classical properties. At the end of each 
section, we will then compare with the exact quantum 
mechanical result. 



III. SEMICLASSICAL THEORY IN THE 
PRESENCE OF A CONTINUOUS SYMMETRY: 
THE ISOTROPIC HARMONIC OSCILLATOR 

As we have seen, the semiclassical density of states is 
intimately related to the classical periodic orbits. Be- 
cause the incommensurate harmonic oscillator exhibits 
only two such orbits, its semiclassical density of states is 
relatively simple to determine. 



As the next step toward a semiclassical understand- 
ing of the commensurate anisotropic harmonic oscillator, 
we consider the more complex case of equal frequencies, 
u)\ = u>2 = UJ- The fundamental periods of oscillations 
along the q±- and c^-axes are equal, T\ = T 2 = T = 2tt/uj 
and from Eq. |T3D we find X 2 (rT 1 ) = Xi(rT 2 ) = 1 for any 
integer r, which is related to the fact that any starting 
point on an energy surface H = E in the four-dimensional 
phase space leads to a periodic orbit. In particular, the 
oscillations along the q^- and g 2 -axes are still periodic and 
should contribute to the semiclassical density of states 
as was the case in Sec. II. However, the denominators 
in Eq. ([8]) reflecting the stability of the orbits vanish. 
The resulting divergences indicate that the evaluation 
of the propagator within the stationary phase approxi- 
mation fails because the stationary points are no longer 
isolated but form a continuous family. A small perturba- 
tion of the initial conditions thus leads to another peri- 
odic orbit close to the original one. It is even possible to 
continuously deform the two Lissajous curves depicted in 
Fig.[IJa) into each other. The existence of a family of pe- 
riodic orbits is a consequence of the presence of a contin- 
uous symmetry that maps a given orbit to a neighboring 
one. Hence, the relation between the classical properties 
and the density of states discussed in the Sec. II has to 
be modified. 

To understand the necessary changes, we consider the 
phase space symmetries of the isotropic harmonic oscil- 
lator and define the quantitie s 15 ' 16 



Ji 



1 / P1P2 
2oj V in 
1 



mcu qiq 2 



•h = ^{qiP2 - qip\) 



7 1 (p\ , m 2 2 p\ m 2 2 



2lo \2m ' 2" 11 2m 2 
for which we find the Poisson brackets 



(16a) 
(16b) 
(16c) 

(17) 



where eijk is the completely antisymmetric tensor. Two 
of the Ji have an obvious physical meaning. J 2 is propor- 
tional to the angular momentum and generates rotations 
around the axis perpendicular to the q\-q 2 plane. It thus 
accounts for the rotational symmetry of the isotropic os- 
cillator in real space. J3 represents the difference of the 
energies of the normal modes. J 2 and J3 are conserved 
quantities as can be verified explicitly by evaluating the 
Poisson bracket with the Hamiltonian of the isotropic 
harmonic oscillator. Because J\ results from the Poisson 
bracket of two conserved quantities, it is a constant of 
the motion as well. We also have 



Jt + Ji 



From the quantum version of Eq. (1171 



[Ji, Jj] = iheijkJk 



(18) 



(19) 



4 



the relation to spin 1/2 is evident. Both the spin 1/2 and 
the isotropic two-dimensional harmonic oscillator exhibit 
a SU(2) symmetry, which indicates how the periodic or- 
bits of the isotropic harmonic oscillator can be classified. 

It is useful to recall the concept of a Bloch sphere. If 
we introduce spherical coordinates, every state of a spin 
1/2 can be expressed in terms of the orthogonal basis 
states |f) and ||) as 

|V) = cos(^/2)|T> + sin(0/2) e ^U), (20) 

where 9 runs from to 7r and ip from to 27r. This rep- 
resentation provides a one-to-one mapping between the 
spin states and the points on a sphere. In particular, the 
cartesian coordinates Xi of a point on the Bloch sphere 
are proportional to the spin expectation values (Si). We 
also note that the spin operators serve as generators of 
rotations on the Bloch sphere so that an arbitrary state 
can be generated by starting from a given initial state, 
for example, |f). 

In close analogy, there is a one-to-one correspondence 
between a point f on a unit sphere and a classical peri- 
odic orbit of the isotropic harmonic oscillator. The or- 
bit is uniquely characterized by the values of the con- 
served quantities in Eq. (|16p which can be expressed by 
J = (E/2<j)r. According to Eq. (I16c|) . the latitude on 
the Bloch sphere determines the distribution of the en- 
ergy between the two oscillation modes along the coordi- 
nate axes. The north and south pole of the Bloch sphere 
correspond to extremal values of J3, that is, to an os- 
cillation along one of the coordinate axes. The equator 
represents all periodic orbits where the energy is equally 
shared between the two modes. The position on the equa- 
tor determines the relative phase of the two oscillatory 
motions. 

In four-dimensional phase space these three conserved 
quantities leave one degree of freedom undetermined, 
which corresponds to the motion along the periodic orbit 
generated by the Hamiltonian. In analogy to spin 1/2, 
all periodic orbits can be generated from a given periodic 
orbit by means of the conserved quantities in Eq. (|16p . 

A generalized trace formula that takes into account 
continuous symmetries was given in Refs. [l7lfl8lfl9l If 
we adapt Eq. (5.1) in Ref. [l8| to the special case of an 
isotropic harmonic oscillator where all periodic orbits are 
members of a single two-parametric family, we find that 
the deviation from the classical density of states of the 
semiclassical (sc) density of states is 

As before, Eq. l[2"T|) is valid only for positive energies. 
It resembles Eq. ((5} in several respects, but there are 
differences due to the existence of a family of periodic 
orbits. We first remark that all members of the family 
share the same action W, period T, and Maslov indices 
<7 r for r repetitions of the periodic orbit. Therefore, apart 



from the factor T — 2ir/uj, which arises from the integral 
along a periodic orbit, an additional factor V\J\ appears 
as a result of the integration over the family. The vol- 
ume occupied by the family is given by the surface of the 
Bloch sphere V = iir. The factor |J| = E/2uj reflects 
our choice for the basis of the Lie algebra characterizing 
the symmetry group A comparison of Eqs. ((SJ and 
(12ip reveals that the factor derived from stability consid- 
erations in Eq. (|15| is missing for the isotropic harmonic 
oscillator. This difference should be expected because ev- 
ery perturbation of a periodic orbit will lead to another 
periodic orbit so that in contrast to the incommensurate 
harmonic oscillator, the question of lateral stability for 
the isotropic harmonic oscillator does not make sense. 

We note that the change in the symmetry properties 
prevents us from deriving the Maslov index from the ex- 
pression Eq. ([9]) by simply setting uj\ = u>2- The result for 
the incommensurate harmonic oscillator has to be modi- 
fied by an additional symmetry-induced contribution^^ 
Act = — 2 so that the Maslov index of the isotropic har- 
monic oscillator becomes 

o r = 4r - 1. (22) 

Substituting the expressions for the period T, the sur- 
face of the Bloch sphere V, and the factor | J\ given be- 
low Eq. (|2"Tj) as well as the action W = 2irE/u> and the 
Maslov index (|22p in the generalized trace formula (|2ip . 
we arrive at the semiclassical correction to the classical 
density of states for the isotropic harmonic oscillator 

2E ^ , 2ttE\ 
E 00 E , (23) 

n— — 00 

In the second step, we have used the Poisson resumma- 
tion formula. After adding the classical density of states 
([3]) we recover for positive energies the exact quantum 
mechanical result as was the case for the incommensu- 
rate harmonic oscillator. 



IV. SEMICLASSICAL THEORY FOR THE 
COMMENSURATE ANISOTROPIC HARMONIC 
OSCILLATOR 

A. From the anisotropic to the isotropic oscillator 

We now turn to the case of different but commen- 
surate frequencies where isolated periodic orbits coexist 
with a family of periodic orbits. A key point for under- 
standing the properties of a commensurate anisotropic 
harmonic oscillator, both with respect to group theory 
and to its semiclassical treatment, is the mapping to the 
isotropic harmonic oscillator. With the spatial symme- 
tries in mind, it might come as a surprise that such a 
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mapping should be possible because an anisotropic har- 
monic oscillator does not obey the rotational symmetry 
characteristic for the isotropic case. In fact, real space 
symmetries are not useful here and we will be required 
to consider symmetries in phase space. 

To motivate the relation between the commensurate 
anisotropic harmonic oscillator and the isotropic har- 
monic oscillator, we consider the energy spectrum of a 
harmonic oscillator with commensurate frequencies sat- 
isfying k\UJ-y — k 2 uJ2 — where k\ and k 2 are assumed to 
be integers without a common divisor. It is instructive to 
express the quantum numbers Vi appearing in the energy 
spectrum ((4]) in the commensurate case as i>i = riiki + A.; , 
where A; = 0, 1, . . . , ki — 1 and rii — 0, 1, 2, . . . ,— so that 
the energy spectrum takes the form 

+ „ + £ + * + * + (24) 

These eigenenergies can be viewed as describing k\k 2 sets 
of spectra of a two-dimensional isotropic harmonic os- 
cillator shifted in energy relative to each other so that 
the isotropic harmonic oscillator appears to be repre- 
sented multiple times in the commensurate anisotropic 
harmonic oscillator. This interpretation motivates us to 
seek a mapping between these two systems and also to 
find their classical counterparts. 

Because we are interested in the phase space struc- 
ture, we briefly recall the results of the Hamilton- Jacobi 
method applied to the harmonic oscillator The action 
Wi is given by Eq. (jTTJ) and the corresponding angle can 
be expressed as 

$i = arctan ( — ; j _ ^5) 

The motion of the harmonic oscillator thus takes place 
on a torus that is parametrized by the angles $i and 
The size of the torus is determined by the actions W\ and 
W2 ■ The solid line in Fig. [Ha) depicts the motion on the 
torus for a commensurate anisotropic harmonic oscillator 
with k\ = 5 and ki = 4. 

According to Eq. (fTTj) , the mapping from a 
commensurate anisotropic harmonic oscillator to an 
isotropic harmonic oscillator can be achieved by the 
transformational^ 

W i ^W i = ^ i (26a) 
^ (26b) 

The definition of the new angles "di ensures that the trans- 
formation is canonical. Remarkably, the different angles 
•di + \i2ir /ki, Aj = 0, 1, . . . ki — 1 are mapped onto the 
same angle fli, because we have to identify di + 2n with $i 
as usual. The transformation (|2"6"|) thus leads to a folding 
of phase space into ki sheets for each mode i. This folding 
is necessary to turn the commensurate anisotropic har- 
monic oscillator into an isotropic harmonic oscillator, be- 
cause it ensures that the original periods Ti = 2i:ki of the 




FIG. 2: The motion of a commensurate anisotropic harmonic 
oscillator with ki — 5 and — 4 is depicted in action-angle 
variables and in spatial coordinates, (a) The full line repre- 
sents a trajectory of the commensurate anisotropic harmonic 
oscillator on a torus in phase space parametrized by action- 
angle variables. The white section indicates one of the fcifc2 
elementary cells delineated by the dashed lines. Each cell 
after identification of its opposite borders corresponds to an 
isotropic harmonic oscillator, (b) The projection of the tra- 
jectory onto real space leads to a Lissajous curve. The small 
circles indicate the positions where the trajectory crosses from 
one elementary cell into another one. The full line represents 
the trajectory in one such cell. 



normal modes are mapped onto equal periods T = 2ir. 
After the canonical transformation (|26p we now have k\ki 
isotropic harmonic oscillators instead of one commensu- 
rate anisotropic harmonic oscillator. Figure [2ja) illus- 
trates the transformation (|2l))) . The white cell covering 
the interval 2ir /ki for angle i?j turns into a square of side 
length 27T by the transformation (|26|) . By identifying op- 
posite sides of the cell, a new torus is formed on which 
the motion of an isotropic harmonic oscillator in phase 
space takes place. Because the original torus contains 
k\kz such cells marked by the dashed lines, Eq. (f2"6"|) pro- 
vides a mapping from a commensurate anisotropic har- 
monic oscillator to k\ki isotropic harmonic oscillators. 



B. Symmetries of the commensurate anisotropic 
harmonic oscillator 

We can now analyze the symmetries present in the 
commensurate anisotropic harmonic oscillator: There ex- 
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ists a discrete symmetry associated with the mapping 
from the commensurate anisotropic harmonic oscillator 
to multiple copies of isotropic harmonic oscillators. In 
addition, according to Sec. IIII1 each isotropic harmonic 
oscillator possesses the continuous symmetry SU(2). The 
presence of the discrete symmetry is a direct consequence 
of the folding of phase space discussed in Sec. IIV Al We 
can introduce complex variables z, = \fWl exp(— to 
write the transformation (12611 as 



N 



\\z,\J ' 



(27) 



which is invariant under Zi t— > &t cxp(2iri / ki) . Hence, 
the symmetry operations consist of rotations in the com- 
plex plane by multiples of 27r/fej. The rotations form a 
symmetry group d = {c gji , cj^, . . . , c^J generated by 
the fundamental rotation c gj i by 2n/ki. The identity is 
given by the group element which leads to a rotation 
by 2tt. Taking together the rotations for both oscilla- 
tor modes, we arrive at the product group C = C\ x C 2 
containing k\k 2 elements and generated by the element 
c g = c g.i x c g,2- For later use we note that the irre- 
ducible representation s 21 ' 22 of C can be labeled by a tu- 
pel (Ai, A2) with Xi = 0,1, ... ,ki — 1. The character 21 
X of the group clement c g in the representation (Ai,A2) 
is given by ^Ai,A 2 (c g ) = exp(27riCA 1 .A 2 ); where we have 
introduced the phases 



r - Al 4- 2 
«i k 2 



(28) 



We emphasize that here we are concerned with a phase 
space symmetry. In real space, a symmetry operation 
would map a segment of a Lissajous curve such as the 
solid curve in Fig. [2Kb) onto one of the other segments 
(delimited in the figure by two small circles) of the same 
Lissajous curve. By starting from a fundamental seg- 
ment, we can arrive at the complete Lissajous curve by 
successively applying the generator c g of the group k\k 2 
times. In the Bloch sphere picture, each segment is rep- 
resented as a point on a Bloch sphere corresponding to 
one isotropic harmonic oscillator. Because the commen- 
surate anisotropic harmonic oscillator has been mapped 
onto kikz isotropic harmonic oscillators, a complete Lis- 
sajous curve corresponds to one point on each of k\k 2 
Bloch spheres. The mapping between two segments of 
a Lissajous curve then implies a jump from one Bloch 
sphere to another one. 



C. Construction of the semiclassical density of 
states 

For the semiclassical treatment of the commensurate 
anisotropic harmonic oscillator, we need to know which 
periodic orbits come in families and which ones are iso- 
lated. We start with a general Lissajous curve where 
both one-dimensional harmonic oscillators along the q\- 



and <72-axes have to return to the starting point of their 
oscillations at the same time. Because their fundamen- 
tal periods are equal to Tj = 2i:ki/uj and k\ and k 2 are 
assumed to have no common divisor, the fundamental 
period for a Lissajous curve is given by T — 2-Kkik 2 /uj. 
A Lissajous curve in configuration space corresponds to 
a periodic orbit in full phase space where Xi(27rfci&2) = 
Y^ 2 (2~Kk\k 2 ) = 1. By varying the energy difference be- 
tween the two oscillators and the relative phase of their 
oscillations, all possible Lissajous curves are obtained, 
including their degeneracies into oscillations along the 
coordinate axes. 

During the time T, k 2 oscillations along the q\ axis 
and ki oscillations along the q 2 axis occur. Hence, all os- 
cillations along the q\ (52) axis with fundamental period 
T\ (T 2 ) that are not repeated k 2 (k\) times or multiples 
thereof are not members of the family of Lissajous curves. 
They cannot be deformed continuously into a Lissajous 
curve nor into each other and lie isolated in phase space. 
Accordingly, for the isolated periodic orbits, one of the 
matrizants Xi or X2 becomes equal to unity but not both. 

These different types of periodic orbits lead to a de- 
composition of the semiclassical density of states of the 
commensurate anisotropic harmonic oscillator into four 
terms 

Psc (E) = p{E) + 5p v (E) + Sp! (E) + 6 P2 (E), (E > 0) 

(29) 

where p is the classical density of states ([3]) . 5pr accounts 
for the contribution of the family of Lissajous curves de- 
noted by r here and in the following, and the last two 
terms arise from the contribution of the isolated periodic 
orbits. The latter can be treated in the same way as the 
two periodic orbits of the incommensurate harmonic os- 
cillator in Sec. [Hi From Eq. ([8]) we find after substituting 
the corresponding period and the Maslov index |9j): 



r—1 



sm 



iy 



sin(7rrfci//c2) 



(30) 



The prime indicates that the sum is restricted to the iso- 
lated periodic orbits, that is, those repetitions r for which 
the denominator does not vanish; 5p 2 (E) is obtained by 
interchanging k\ and k 2 in Eq. (|30p . 

To determine the contribution Sp^ of the family of Lis- 
sajous curves to the density of states, we make use of the 
discrete symmetry presented in Sec. IIVB1 The mapping 
to the isotropic harmonic oscillator will allow us to take 
advantage of the results of Sec. Mil 

The presence of a symmetry imposes a structure on the 
dynamical equations that lets us consider the kinematics 
in a reduced space. In a quantum mechanical picture, 
the full Hilbert space of the problem can be decomposed 
into invariant subspaces corresponding to the irreducible 
representations of the underlying symmetry group— As 
a consequence, the solution of the Schrodinger equation 
that is restricted to these subspaces yields symmetry- 
reduced spectra. 
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A semiclassical treatment of this idea was used in 
Refs.l23T-[25lfor discrete symmetries in configuration space 
and in Ref. [26| for continuous symmetries. Although the 
commensurate anisotropic harmonic oscillator exhibits a 
discrete as well as a continuous symmetry, we use symme- 
try reduction only for the discrete group C. After we have 
determined the symmetry reduced densities of states, we 
will sum over all irreducible representations. It can be 
generally shown that the full spectrum is constructed in 
this way32£>2£ 

An orthogonal projection corresponding to the irre- 
ducible representation (Ai, A2) can be achieved using the 
projection operator— 



hup(E) 



MM) 



(31) 



where U(c) is the unitary operator describing the action 
of the group element c in Hilbert space. g?Ai,A 2 is given by 
the dimension of the irreducible representation which in 
our case is always equal to one^ By this projection the 
Green function G(E), its trace and, in view of Eq. (fTUJ), 
the density of states arc modified accordingly. 20 The re- 
duced densities of states corresponding to the irreducible 
representation (Ai,A2) are found to be 



5p ( r XlM \E) 
2E 



(32) 



h 2 u 



2nr 



The phases Cai,a 2 are defined in Eq. (|28[) and arise from 
the character X present in Eq. (f3Tj) . As we will see in 
the discussion of Fig. (3[ this phase is responsible for a 
relative shift of the reduced energy spectra. 

The expression (|32l) can also be interpreted within a 
classical picture, where the available phase space is re- 
duced to a primitive cell, for example, the white area de- 
picted in Fig.[2)Ja). The motion in full phase space may be 
recovered using symmetry operations. Correspondingly, 
a Lissajous curve of the commensurate anisotropic har- 
monic oscillator is reduced to a fundamental segment, for 
example, the full line in Fig.^b). The endpoints of such 
a segment are identified, thus yielding a periodic orbit of 
the isotropic harmonic oscillator whose continuous SU(2) 
symmetry still remains. We therefore can directly apply 
the results of Sec. IIIII It turns out that the contribution 
to the semiclassical density of states of the rth repetition 
of such an orbit acquires an additional phase given by the 
character of the group element Cg. This group operation 
is needed to close the corresponding partial orbit in full 
phase space and, within a classical picture, can be viewed 
as the origin of the shift of the reduced energy spectra 
mentioned previously. 

It remains to determine the Maslov index 5y of an orbit 
segment in Eq. ([32)) . Let us consider as a representative 
of the whole family the oscillation along the qi axis. We 
recall that only k 2 repetitions of the fundamental peri- 
odic orbit form a member of the family. We also have 




FIG. 3: Density of states of the commensurate anisotropic 
harmonic oscillator with frequency ratio 3:2. (a) Contribution 
of the classical density of states and the family of periodic 
orbits. The six different symmetry classes (Ai, A2) that repeat 
themselves periodically are indicated, (b) The contribution 
arising from the isolated periodic orbits, (c) The total density 
of states of the commensurate anisotropic harmonic oscillator. 



to keep in mind that we are now working in reduced 
dynamics where the k\k 2 th part of a periodic orbit of a 
commensurate anisotropic harmonic oscillator appears as 
a periodic orbit of an isotropic harmonic oscillator. We 
therefore have to replace r in Eq. ([9|) by rk 2 /k\k 2 . If 
we add the correction Act = —2 as in the derivation of 
Eq. (fJU) for the isotropic harmonic oscillator, we obtain 



5y = 2r 



fci + k 2 
k 1 k 2 



1. 



(33) 



With Eqs. ©, ([3D]), (J5J), and (J35J we have all the in- 
formation needed to construct the full density of states 
([29]) . We proceed in two steps as illustrated in Fig. [3] for 
the special case of k\ = 3 and k 2 — 2. First, we note 
that the symmetry reduction performed for the family of 
periodic orbits can also be applied to the classical density 
of states p(E). As we have seen, the reduced dynamics 
takes place in only the kik 2 th part of the original phase 
space volume. Correspondingly, the classical density of 
states contributes to each reduced density of states only 
p{E) jk\k 2 . We add this contribution to each of the re- 
duced densities of states arising from the family of peri- 
odic orbits and arrive at k\k 2 series of delta peaks with 
a height increasing linearly with energy. The density of 
states resulting from these two contributions is shown in 
Fig. [3]^a) where k\k 2 = 6 irreducible representations lead 
to a corresponding number of series of delta peaks shifted 
relative to each other in energy. 

This density of states has to be corrected by the contri- 
bution Spi(E) + 8p 2 (E) arising from the isolated orbits 
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(see Fig. [31(b)). The sum of all contributions is shown 
in Fig. [3Jc) , where we can check that the eigenenergies 
with their degeneracies are correctly generated. We re- 
mark that the first two peaks of Fig. Ela) are canceled 
by the first two peaks in Fig. [3fb), so that the ground 
state energy 5fiw/12 is associated with the irreducible 
representation (0,0). 

Although this example indicates that the semiclassical 
treatment reproduces the spectrum of the commensurate 
anisotropic harmonic oscillator, it is useful to explicitly 
perform the sum over the reduced densities of states ([3"2"1) 

fci-l k 2 -l 

Spr(E) = £ E ^r A " A2) (£)- (34) 

Ai=0 A 2 =0 

By interchanging the double sum with the sum over all 
repetitions r in Eq. (|32|). it is seen that Eq. ([34]) vanishes 
except when r is a multiple of k\k2- Therefore only such 
orbits of the reduced dynamics actually contribute, which 
are repeated k\k2 times. As expected, these orbits are 
exactly the full Lissajous curves. 

We may now replace the summation index r in Eq. (|32[) 
by r = k\k2k where k = 1, 2, . . ., and arrive at 

Spr(E) (35) 
fc=i 

Adding the classical density of states and the contribu- 
tions of the isolated orbits, Eq. ([29]) yields the exact den- 
sity of states of the commensurate anisotropic harmonic 
oscillator £ This result can be checked by summing the 
exact partition function and taking the inverse Laplace 
transform as we did for the incommensurate harmonic 
oscillator in Sec. [H] 

V. CONCLUSIONS 

The harmonic oscillator is a unique system with funda- 
mental applications in many areas in physics. It admits 
exact solutions that are often helpful in the illustration 
of theoretical approaches and methods. As we have seen, 
it is particularly interesting to consider a semiclassical 
treatment where a relation between classical and quan- 
tum properties can be established. Beyond the interest- 
ing different scenarios for periodic orbits provided by the 



two-dimensional harmonic oscillator, this system beauti- 
fully illustrates the relation between classical phase space 
symmetries and quantum degeneracies. 



VI. SUGGESTED PROBLEMS 

Problem 1. The conserved quantities Jk of the isotropic 
harmonic oscillator defined in Eq. (|16p can be interpreted 
as generators of infinitesimal transformations. For exam- 
ple, J2 generates rotations around the axis perpendicular 
to the Xi-%2 plane. To better understand the significance 
of Ji, consider its action on the energies Hi and H2 of 
the two modes, that is, evaluate the Poisson brackets 
{Hi, Ji}. How does J\ modify an orbit? Hint: The ratio 
of H\ and H2 is related to the eccentricity of the orbit, 
see Ref. [l6| for a discussion. 

It is also instructive to consider the Jk as quantum- 
mechanical operators. Express them in terms of creation 
and annihilation operators b\ and 6j, respectively, for the 
two modes. Determine the action of J\ ±i J2 on an eigen- 
state \v\,v 2 ) with Hi\v\,vi) = fiw(vi + 1/2) j 1^1 , v 2 ). 

Problem 2. For the one-dimensional harmonic oscilla- 
tor a transformation of the density of states analogous to 
but simpler than the transition from Eq. ^ to Eqs. |7|)- 
§§§ can be performed. From the partition function Zp of 
the one-dimensional harmonic oscillator, determine the 
density of states by an inverse Laplace transformation 

p(E) = — d0Zf,exp{l3E), (36) 

2mJ t _ loo 

where e > 0. Use residue calculus and close the contour 
in the left complex half plane. The resultant expression 
for the density of states should contain the classical con- 
tribution p = 1/fvu) and an oscillating part 

r=l 

How can this result be interpreted in terms of the classical 
properties of the one-dimensional oscillator? Discuss the 
differences between the one- and two-dimensional case. 
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